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ABSTRACT: A Monte Carlo event generator is constructed for a two-Higgs-doublet model 
with maximal CP symmetry, the MCPM. The model contains five physical Higgs bosons; the 
p' , behaving similarly to the standard-model Higgs boson, two extra neutral bosons h! and 
h", and a charged pair . The special feature of the MCPM is that, concerning the Yukawa 
couplings, the bosons h! , h" and couple directly only to the second generation fermions 
but with strengths given by the third-generation-fermion masses. Our event generator 
allows the simulation of the Drell-Yan-type production processes of h', h" and in proton- 
proton collisions at LHC energies. Also the subsequent leptonic decays of these bosons 
into the // + //~, p + v^ and [i~Vn channels are studied as well as the dominant background 
processes. We estimate the integrated luminosities needed in pp collisions at center-of-mass 
energies of 8 TeV and 14 TeV for significant observations of the Higgs bosons h', h" and 
in these muonic channels. 
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1 Introduction 

The experimental investigation of the Higgs sector of particle physics is one of the main 
aims of the LHC experiments. In the standard model (SM) there is only one physical Higgs 
boson. But more complicated Higgs sectors are by no means excluded. In this article 
we present a Monte Carlo event generator for a particular two-Higgs-doublet model. The 
construction of such an event generator is necessary for allowing realistic comparisons of 
theory and experiment to be done for this model. 

Our paper is organised as follows. In section 2 we recall the main construction principles 
and properties of the "maximally CP-symmetric model" (MCPM) which will be studied. In 
section 3 we present predictions for the LHC experiments. In section 4 the implementation 
of the MCPM into the Monte Carlo event-generation package MadGraph and its validation 
are summarized. Section 5 deals with a brief analysis of MCPM signatures using the new 
Monte Carlo tools. We draw our conclusions in section 6. 

2 The Maximally CP-Symmetric Model 

Extending the Standard Model Higgs sector to two Higgs doublets, 

-ft)' "-{%)• 
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gives the two-Higgs-doublet model (THDM). There, the potential may contain many more 
terms than in the SM; see e.g. [1, 2]. The most general THDM Higgs potential can be 
written as follows [3] 

V = ml x ((p\<pi) + m| 2 (4<p 2 ) - m 2 12 ((f\(p 2 ) - {m\ 2 )* {^i) 

+ ^l(yi^l) 2 + ^ 2 (4^2) 2 + \z{<p\yi){<p\<P2) 

+ [Mvlw) + W<P2<Pi)K<Pi<Pi) + iM^m) + a?(4vi)](4v2) . ( 2 - 2 ) 

with rriii, 77J 22 , Ai 2,3,4 real and m\ 2 , ^5,6,7 complex. Many properties of THDMs turn out 
to have a simple geometric meaning if we introduce gauge invariant bilinears [4, 5], 



K = (p\ipi + ip 2 (p2, K 



\K 3/ 



( <Pi<P2 + <fl<Pl \ 



up\yi - if{(f2 ■ (2.3) 



In terms of these bilinears Kq, K, the Higgs potential (2.2) reads 

V = £ K + £ T K + m K$ + 2K oV T K + K T EK (2.4) 

with parameters £q, 770O; 3-component vectors £, r\ and a 3 x 3 matrix E = E^, all real. 
The standard CP transformation of the Higgs-doublet fields is defined by 

<Pi(x) ->■ <p*( x '), i = l,2 x' = (x°,-x). (2.5) 

In terms of the bilinears, this standard CP transformation is [6, 7] 

K (x) -»• Ko(x') , K(x) -»• R 2 K(x') (2.6) 

where R 2 = diag(l, — 1, 1), corresponding in K space to a reflection on the 1-3 plane. 
Generalised CP transformations (GCPs) are defined by [8-10] 

tpiW-tUijipfa), i,j = 1,2, (2.7) 

with U being an arbitrary unitary 2x2 matrix. In terms of the bilinears this reads [7] 

K (x) -> K (x'), K{x) -> R K(x') (2.8) 

with an improper rotation matrix R. 

Requiring R? = I3 leads to two types of GCPs. In K space: 

(i) R = — 13, point reflection, (2-9) 

(ii) R = R T R 2 R, reflection on a plane (R G 5*0(3)). (2.10) 

For a review of the bilinear formalism and its relation to the conventional field approach 
as well as the generalised CP transformations we refer to [11]. There, also an extensive list 
of references is given. 
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While the CP transformations of type (ii) in (2.10) are equivalent to the standard CP 
transformation (2.5), respectively (2.6), up to a basis change, the point reflection trans- 
formation of type (i) is quite different and turns out to have very interesting properties. 
Motivated by this geometric picture of generalised CP transformations, the most general 
THDM invariant under the point reflection (i) has been studied in [12-15]. Invariance of 
the THDM potential (2.4) under the GCP transformation (2.8) with R = -1 3 from (2.9) 
clearly requires £ = T] = 0, leading to 

^MCPM = Co^o + r?oo^o+ KT ^ K - (2-H) 
Without loss of generality the 3x3 matrix E can be chosen to be diagonal 

E = diag(/xi,/i 2 ,^3) (2.12) 

with the ordering 

Ml > M2 > fJ-3- (2.13) 
The conditions to obtain a physically acceptable theory are spelled out in [12, 13] and read 

f?oo > , 
Ma + %0>0, for a = 1,2, 3, 

(2.14 

£o<0, 
/*3 < 0. 

The potential (2.11) of this model is, besides the point reflection symmetry of type (i), 
invariant under three GCPs of type (ii); see (2.9), (2.10). Requiring also the Yukawa 
couplings to respect these four GCPs, it was found in [12] that at least two fermion families 
are necessary in order to have non- vanishing fermion masses. That is, a reason for family 
replication was given. For two fermion families only three options for the Yukawa-coupling 
structure were found: 

(a) one massive and one massless family, 

(b) two mass-degenerate families, 

(c) large flavour changing neutral currents. 

The maximally CP-symmetric model (MCPM) was constructed as follows in [12, 13]. The 
option (a) above was chosen for the Yukawa couplings with the third family (t, b, r) as 
massive one, the second family (c,s,fi) as massless one. The first family (u,d,e) was 
added uncoupled to the Higgs fields. In this way the highly symmetric MCPM gives a 
first approximation to the fundamental mass spectrum as observed in nature; see section 5 
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of [12]. The Yukawa couplings of the MCPM read 

{««> 

-S{^),i(x)(^) 

- S - R (,)^)|^jjJ+k. (2.15) 

where c[g , c^ 1 ] and cjjg are real positive constants, determined by the fermion masses as 
discussed below. 

In the unitary gauge electroweak symmetry breaking (EWSB) gives 

Vl{x) = 7=2 U ' ^ = (^m+wix))) (2 - 16) 

with the standard vacuum expectation value 



v = J — « 246 GeV. (2.17) 

V ^oo + to 

The physical Higgs-boson fields are pf, ti, h", H+, and H~ = (H + ) j . Their masses - at 
tree level - read [12, 13] 

m 2 p , = 2(-£ ) , 

m l =2^o(Mi-M3), 

""V = 2 ?W2 ~ to) , 

m 2 H± = 2v 2 (-n 3 ) 

where the ordering of the masses 

m 2 h , > m 2 h „ (2.19) 

is predicted by the theory; see (2.13). Requiring now that the neutral Higgs bosons b! and 
h" are neither massless nor mass degenerate leads to the condition 

Ml > to > to, (2-20) 
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replacing the weaker one in (2.13). The masses of r, t, b are given by 



m T = c l3 



™t = c^^=, (2.21) 



(l) 

v/2' 

(1) 

(l) v 
m b = c d3^ ■ 

This fixes the constants of the Yukawa couplings in (2.15) in terms of the third- family 
fermion masses. 

The complete Lagrangian after EWSB is given in appendix A of [13] and is implemented 
as Monte Carlo program for the Feynman rules in appendix A of the present paper. 

In total, there are 11 parameters in the MCPM Lagrangian. Before EWSB the natural 
parameters are the gauge couplings g, g' , g s , the Higgs potential parameters p\, p2, P3, 
£o, rjoo, and the Yukawa coefficients cj^ , c^\, cQ. It is convenient to replace the above 
11 parameters by another set of 11 independent parameters containing more directly and 
well measurable quantities. These are the electromagnetic and strong coupling constants 
a and a s , respectively, the Fermi constant Gf, the Z-boson mass rriz, the fermion masses 
mt, mb, m T , and the Higgs boson masses m p i, m^/, my/, m H ±. The relation of these 11 
parameters to the original ones is given in appendix A. The masses of the physical Higgs 
bosons are constrained by (2.19) and from the experimental results on the so called oblique 
parameters in electroweak precision tests; see [16]. 

Let us summarize the essential properties of the MCPM: 

• There are 5 physical Higgs particles, two CP even ones p' , b! , one CP odd one h" , 
and a charged Higgs-boson pair . 

• The p' boson's fermionic couplings are exclusively to the third (r, t, b) family. The p' 
behaves similarly to the SM Higgs boson. 

• The Higgs-boson-fermion couplings of h', h", are exclusively to the second (p, c, s) fam- 
ily with strengths proportional to the masses of the third generation fermions. 

• The first (e, u, d) family is uncoupled to the Higgs bosons. 
For further details we refer to [12-16]. 



3 Predictions for Hadron Colliders 

Since the Yukawa couplings of the h', h", Higgs bosons to the second fermion family 
are proportional to the third-fermion-family masses we have for them large cross sections 
for Drell-Yan type production. For the same reason we have large decay rates of these 
Higgs bosons to the second generation fermions. In figure 1 we show the diagrams for 
these production and decay reactions for the h', h", and H~ bosons in pp collisions. The 
corresponding diagram for H + production and decay is similar to the right one of figure 1 
with the replacements H~ — > H + , s — > c, c — > s, p~ — > and z7„ —> p + . In [13, 15] the cross 
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Figure 1: Feynman diagrams for the dominant Drell-Yan type production processes for 
the Higgs bosons ft, h" and H~ in the MCPM (q = c,s). The dominant hadronic and 
leptonic decays are also indicated. 
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Figure 2: Total cross section of Drell- 
Yan type Higgs boson production at the 
LHC for c. m. energies of 8 TeV and 
14 TeV, respectively. 
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Figure 3: Branching ratios of the CP 
odd h" Higgs boson, where a mass of 
mjj± = 200 GeV is assumed. Figure 
taken from [13]. 



sections were computed for Drell-Yan Higgs-boson production at the LHC for center-of- 
mass energies of 7 TeV and 14 TeV, respectively. In [14] radiative effects were considered. 
Here, we add the cross sections for a center-of-mass energy of 8 TeV which is currently 
available at the LHC. The corresponding total cross sections for the Drell-Yan production 
of the ft ', h", bosons are shown in figure 2. In figure 3 we also recall the branching ratios 
of the h" boson decays. In the following we shall develop all the necessary theoretical tools 
for detailed comparisons of the MCPM predictions with experimental data at the LHC. In 
particular, a Monte Carlo event generator for the MCPM shall be constructed. 

4 Implementation of Monte Carlo Event Generation 

We have implemented the MCPM into the Monte Carlo event-generation package Mad- 
Graph 5 / MadEvent and validated this implementation with different methods. 

MadGraph 5 [17] is the latest version of the matrix-element generator MadGraph. It is 
bundled with the event-generation package MadEvent. For simplicity, both will be referred 
to as MadGraph in the following. Given an arbitrary process for any implemented model, 
the program lists all contributing tree-level diagrams. For each of these diagrams, events 
are randomly produced and the cross section or decay width is calculated. The output of 
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such a simulation is a list of events in the Les Houches Event File (LHEF) standard [18]. 
For each event it includes the particles involved and their four-momenta. This can be read 
by programs such as PYTHIA [19], with which hadronisation effects can be simulated. 

MadGraph is a flexible framework and allows the implementation of new models through 
the Universal FeynRules Output (UFO) format [20]. Its major downside is the limitation 
to tree-level processes. However, the most interesting and relevant processes in the MCPM 
do not involve any loops [13], so this drawback is acceptable for our present study. 

The implementation includes all particles and couplings of the MCPM. This not only 
allows the analysis of signal processes involving Higgs bosons, but also the simulation of 
any tree- level background process. The quantities listed at the end of section 2, namely 

q, Gf, ct s , mz, tth, rnb, m T , m p >, m h >, mh", and m H ± , (4.1) 

are used as independent parameters. 

In addition, more parameters can be used for the study of SM background processes. 
At tree level the MCPM predicts massless fermions for the first and second generation and 
a unit CKM matrix, Vqkm = 13- Of course, this is not what we observe in nature, but it 
may represent a first approximation. Indeed, the ratios of the masses of first and second 
generation fermions to the corresponding third generation ones are small and Vqk M is close 
to the unit matrix; see (125) and (126) of [12]. It may thus be that the GCPs of the MCPM 
play the role of approximate symmetries. Thus, it should be sensible to analyse a model 
that combines the scalar sector of the MCPM with fermions of non-zero mass and a CKM 
matrix Vqkm 7^ I3; as done already in [13]. 

For this reason our MadGraph implementation of the MCPM lets the user set fermion 
masses for the first two families as well. This affects phase space calculations, but the 
interactions remain unchanged. In a similar way the CKM matrix can be set to an arbi- 
trary matrix, which changes SM processes, but not the Higgs boson vertices. With these 
parameters the user can decide whether to analyse the strict MCPM or a hybrid theory 
with fermion masses and a SM-like CKM matrix. 

The MadGraph implementation of the MCPM was validated in two different ways. 
First, total cross sections and decay widths were checked for a number of different pro- 
cesses and parameter sets. These include the production of the MCPM Higgs bosons via 
quark-antiquark fusion as well as their fermionic decays. In all processes and parameter 
configurations good agreement between Monte Carlo results and theoretical expectation 
was found. 

Second, angular dependencies and invariant mass distributions were analysed for dif- 
ferent processes including the Drell-Yan-type production of Higgs bosons followed by their 
decay into fermion pairs. Again, the event shapes agree well with the expected distributions. 
For more details see [21]. 

5 Analysis of MCPM Signatures 

It is now straightforward to ask whether the MCPM can be discovered or excluded at the 
LHC. A full answer requires a thorough analysis including the simulation of fragmentation 
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Parameter set 


m p > 


m h > 


r h > 


m h " 




m H ± 




A 


125 


200 


12.08 


150 


9.06 


150 


9.07 


B 


125 


400 


24.16 


300 


18.12 


300 


18.14 



Table 1: Higgs-boson mass parameters used for the comparison of MCPM signatures with 
SM background. For h', h" and the widths as calculated from table 3 of [13] are 
also given. For the mass parameters chosen these widths are practically equal to the total 
widths. All numbers are in units of GeV. 

and detector behaviour, which goes beyond the scope of this publication. The hadronic 
decay modes will thus not be considered here. A data analysis would neccessiate the 
tagging of charm jets which is experimentally challenging in the presence of huge QCD 
backgrounds. A study using a parton-level MC generator would not yield sensible results. 

However, there is no hadronisation for leptons. Muons, which play an important role in 
the MCPM can be reconstructed relatively precisely in experiments. Hence, for the muonic 
decay channels the hard-process events generated by MadGraph are worth a look even with- 
out fragmentation and detector simulation. As a first application we use the new MadGraph 
implementation to compare the muonic MCPM signatures to the SM background. 

Two such channels are analysed. The first one consists of the production of h! or 
h" via quark-antiquark fusion with subsequent decay into a /U~/i + pair. The dominant SM 
background is given by the 7* and Z Drell-Yan processes. As second channel we analyse the 
production of a H~ boson decaying into a liTv^ final state. Here the dominant background 
stems from the production of W~ bosons. The corresponding Feynman diagrams can 
be found in figure 1. Additional backgrounds from e.g. top production have not been 
considered. Experimental studies show that they are below 15 % for most of the tested 
phase space [22-25]. The production and decay of H + bosons is analogous to H~ bosons 
and could be analysed in the same way. 

Signal and background processes are simulated separately. Due to the different helicity 
structure of the lepton pairs from gauge-boson decays on the one hand and from Higgs-boson 
decays on the other hand there is, nelecting the lepton masses, no interference of signal and 
background. The couplings of [i + [x~ to 7* and Z are chirality conserving. Thus, in the 
limit m^t = only the helicity combinations (/i_|_,/t_) = (5,-3) and (—3, 5) occur for the 

pair. The h! and h" bosons, however, couple in a chirality-changing way to 
Hence, the helicity combinations of fJ, + /J,~ can only be (/&+,/&_) = (g, 5) anc ^ ( — \->~ \)- 
Similar considerations apply to W~ — > H~Vn versus H~ — > V^v^ and W + — > fJ^Vn versus 

Two different parameter sets for the Higgs-boson masses are used, see table 1. For all 
other particles, PDG recommendations [26] are chosen. For each channel and parameter 
configuration, at least 200,000 events (1,000,000 events) are generated for the signal (back- 
ground) processes. All simulations are performed for the LHC design energy ^fs = 14 TeV 
as well as its current center-of-mass energy ^/s = 8 TeV. 

The analysis is first performed without applying any phase-space cuts. Of course, 
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there is no detector that is able to detect and reconstruct every muon. For instance, in the 
ATLAS detector at CERN the muon trigger chambers only cover the pseudorapidity region 
\rj\ < 2.4; see [27]. In addition, selection cuts on the transverse momentum of the muons 
are applied by the experiments in order to guarantee efficient triggering and background 
suppression; a typical cut value is px > 25 GeV. Therefore, we also perform the analysis 
requiring 

M < 2.4, (5.1a) 
p T > 25 GeV (5.1b) 

for each muon. In the /x - ^ channel, 

#p iss > 25 GeV (5.1c) 

is also required, where EJp lss is the missing transverse energy due to the non-detection of 
the neutrino. 

For the channel, the resulting distribution of the invariant mass m(/j,~ , of the 

muon pair is shown in figure 4 for a center-of-mass energy of 1 /s = 14 TeV and in figure 5 
for y/s = 8 TeV. Here and in the following the bin size in mass is chosen as Am = 5 GeV. 
Note that the resonance peaks from the Higgs bosons are enhanced by a factor 10. Thus, 
these Higgs-boson resonances are tiny in comparison to the background. The application 
of selection cuts and a higher center-of-mass energy improve the situation slightly. 

For the fj,~ Vn channel, the event distribution in the transverse mass tut of the lepton 
pair is given in figures 6 and 7 for ^/s = 14 TeV and 8 TeV, respectively. Again, the Higgs- 
boson resonances are much smaller than the background, even using appropriate selection 
cuts and a center-of-mass energy of 14 TeV. 

Using the distributions after the application of selection cuts we calculate how much 
statistics is needed so that the MCPM resonance peaks are locally significant compared to 
the background, i. e. larger than na, where a is the statistical uncertainty of the background. 
Conventionally, an excess with n = 5 is considered a discovery, while n = 2 is often used 
for exclusion limits. 

The detector resolution is a limiting factor in the search for resonances. Depending on 
the channel, the search strategy, the energy scale, the pseudorapidity region and a number 
of other parameters, the ATLAS collaboration quotes e.g. a design value for the dimuon 
invariant-mass resolution between 2 and 12 GeV; for the transverse-mass resolution the 
quoted values are slightly worse [27]. Based on these figures and on the widths V of the 
Higgs bosons as given in table 1, two to five bins were taken into account when calculating 
the integrated luminosities needed for a local 2a- and 5<7-significance, depending on the 
mass scale. 

In table 2 we give the integrated luminosities needed for local 2a- or 5<7-significance 
of the peaks due to the Higgs bosons h', h" and H~ of the MCPM, as shown in figures 4b, 
5b, 6b and 7b. It should be stressed that these results are very rough estimates. On the 
one hand, only the local significance in one channel for the selection cuts given in (5.1) is 
described. The sensitivity of the search might be improved by selection-cut optimisation 
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(b) 

Figure 4: MCPM events compared to the SM background in the channel as a 

function of the invariant mass of the muon pair for pp collisions at a center-of-mass energy 
of yfs = 14 TeV. Two different parameter sets A and B are used for the MCPM, see table 1. 
(a) No selection cuts applied; (b) selection cuts on rj and p? of the muons applied. 



for the different channels. On the other hand, only hard tree-level processes are taken into 
account and the behaviour of the detector was not simulated. 

Keeping this caveat in mind, these numbers may nevertheless hint at the order of 
magnitude of statistics needed for MCPM signatures to become visible. 
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Figure 5: Same as figure 4, but with a center-of-mass energy of yfs = 8 TeV. (a) No 
selection cuts applied; (b) selection cuts on i] and px of the muons applied. 



6 Conclusions 

In this paper we presented the implementation of a two-Higgs-doublet model with maximal 
CP symmetry, the MCPM, into the Monte Carlo event-generation package MadGraph. 

The MCPM is an extension of the Standard Model (SM) that is based on the require- 
ment of invariance under certain generalised CP transformations. It features a fermion 
structure in which only the third-generation fermions are massive. Of course this does 
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Figure 6: MCPM events compared to the SM background in the /x~z^ channel as a 
function of the transverse mass of the final state lepton pair for pp collisions at a center-of- 
mass energy of y/s = 14 TeV. Again the parameter sets A and B are used for the MCPM, 
see table 1. (a) No selection cuts applied; (b) selection cuts on r\ and pt of the muon as 
well as EJp ms applied. 



not describe nature precisely, but the model gives a first approximation of what has been 
observed. The theory predicts five physical Higgs bosons. Concerning their Yukawa cou- 
plings, four of these bosons couple only to the second-generation fermions, but with coupling 
strengths given by the third-generation fermion masses. At colliders such as the LHC, these 
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Figure 7: Same as figure 6, but with a center-of-mass energy of y/s = 8 TeV. (a) No 
selection cuts applied; (b) selection cuts on rj and px of the muon as well as E™ ss applied. 



four Higgs bosons will be produced mainly via quark-antiquark fusion. Their /j> + fj>~, fj vv^ 
and [i~v^ decay modes seem most promising for a discovery. For a proper search for MCPM 
signatures a Monte Carlo simulation is needed. 

We implemented the MCPM into MadGraph, allowing the calculation of cross sections 
and the random generation of events for arbitrary tree- level processes. The output is in 
the LHEF format and can be used for further analysis, for instance using PYTHIA and 
GEANT. Therefore, this MadGraph implementation can be used as the starting point for a 
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Necessary int. luminosity f Ldt [fb x ] 
Parameter set Process ran [GeV J - 

2d significance 5cr significance 



A 


ti' - 




150 


5.2 


33 


A 


ti - 


► n + n~ 


200 


17 


100 


B 


ti' - 


+ //+//" 


300 


99 


620 


B 


ti - 




400 


490 


3100 


A 


H~ 


-> 


150 


9.3 


58 


B 


H~ 




300 


140 


890 



(a) V« = 8 TeV 



Necessary int. luminosity f Ldt [fb x ] 
Parameter set Process ran [GeVJ - 

2d significance 5cr significance 



A 


ti' - 




150 


1.3 


8.3 


A 


ti - 


> n+fjT 


200 


3.4 


21 


B 


ti' - 




300 


17 


104 


B 


ti - 




400 


64 


400 


A 


H- 


-> 


150 


3.1 


20 


B 


H- 




300 


30 


190 



(b) ^fs = 14 TeV 

Table 2: Integrated luminosities that are needed for local significances of the Higgs-boson 
resonances in the invariant mass distribution and transverse mass distribution of the muonic 
final states /U + /i~ and fJ>~Uu, respectively. The figures for H + production and its decay to 
H^Vn, are similar to those for H~ . 



full and thorough Monte Carlo simulation of the MCPM. It was validated successfully with 
different methods. The event generator is available from [28]. 

The implementation was then used to compare the MCPM signatures to the SM back- 
ground at LHC energies. This analysis was only done for the muonic channels and restricted 
to hard processes, so the results are only a rough approximation. For a Higgs boson h" (H~) 
of 300 GeV, an integrated luminosity of approximately 100 fb -1 (140 fb -1 ) at a center-of- 
mass energy y/s = 8 TeV has to be collected so that a local 2ct significance becomes possible. 
On the other hand, for h" (H~) of mass 150 GeV the necessary integrated luminosities at 
y/s = 8 TeV for local 2<r significance are only 5.2 fb -1 (9.3 fb" 1 ). At present (August 2012), 
the ATLAS and CMS experiments have both collected data representing roughly 13 fb -1 
of integrated luminosity. This number is expected to rise up to 20 fb" 1 in the current data 
taking period. Thus, it seems that the exclusion of a part of the Higgs-boson mass range 
or the discovery of small local excesses hinting at the MCPM might be within reach in the 
near future. 
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Finally, we note that the recent announcement [29, 30] of a signal for a boson of mass 
125 GeV is - so far - not only compatible with the SM Higgs boson, but also with the 
SM-like Higgs boson p' of the MCPM. Only detailed comparisons of the decay channels of 
the discovered boson with theoretical predictions will allow to draw further conclusions. We 
also note that the bounds on the masses of the h', h" and bosons given in [16] depend 
on the mass of the p' . Identifying the new boson with the p' of the MCPM, these bounds 
can be sharpened. Only figure 2 of [16] remains relevant where m p i = 125 GeV was chosen. 
Our parameter sets A and B of table 1 also correspond to this choice of the p' mass. 
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A Appendix 

1.1 Feynman Rules of the MCPM 

In appendix A of [13] the MCPM Lagrangian after electroweak symmetry breaking is pre- 
sented and the Feynman rules for several vertices are derived. The most important vertices 
and the corresponding expressions are given in table 3 and table 4 below. 

1.2 Using the Monte Carlo Generator 

The MCPM model in MadGraph is used in the same way as every MadGraph model, see 
[17] for an overview. First it has to be loaded. After starting the MadGraph binary in a 
shell this can be done with the command 

import model MCPM -modelname 

where the option -modelname is needed to ensure the correct particle names. Then processes 
are defined and the event generation is started: 

generate p p > hi > m- m+ 
add process p p > h2 > m- m+ 
output -f 
launch 

The "generate" and "add process" commands are used to specify the processes Mad- 
Graph has to evaluate, while "output" and "launch" start the event generation. The 
particle names used in the implementation are given in the following section. Before events 
are produced, MadGraph asks the user whether the standard parameters are to be used or 
if he wants to modify them. 

During the simulation, an in-browser status page keeps the user informed. In the 
end, MadGraph gives out the contributing Feynman diagrams, the cross-sections or decay 
widths, the corresponding uncertainty and a compressed file containing the generated events 
in LHEF [18] format. This can be used for further analysis, for instance with Mad Analysis 
or PYTHIA [19]. 

1.3 MCPM Parameter Relations 

Here we give the relation - at tree level - between the original parameters of the MCPM 
Lagrangian and the parameters of (4.1) used in the event generator. With vq given in terms 



-17- 



of the original parameters in (2.17) we have 




(A.l 



(A.2 



(A.3 



(A.4 



(A.5 



(A.6 



(A.7 



m 2 p , = 2v 2 (t] 00 + M3) , 
™ 2 h> = 2^oOi-P3), 



(A.8 
(A.9 
(A.10 
(A.ll 



m l" = 2^0(^2-/^3) 
m 2 H ± = 2v 2 (-/i 3 ) . 



1.4 List of Particles and Parameters in the MadGraph Implementation of the 



A list of MCPM particles and parameters with their MadGraph names are presented in 
table 5 and 6, respectively. 

The parameters may be a bit confusing, because the c, e, fi can be set to be massive, 
while the u, d and s are always massless. This does not have anything to do with the MCPM. 
As explained in section 3, some of the consequences of strict symmetry under generalised CP 
transformations have been dropped in this implementation, so it is generally possible to set 
masses for the first- and second-generation fermions. However, MadGraph treats the very 
light fermions as massless (even in the SM, where the theory certainly predicts something 
else). This does not change the physics at collider experiments in the TeV range, but it 
saves some computation time. Therefore the list of parameters includes the masses of those 
fermions of the first two generations where the masses are relevant for the calculations. 

MadGraph also allows to set decay widths for all particles, which are named similarly 
to the mass parameters. This is important for further analysis with programs such as 
PYTHIA; these parameters do not influence the generation of events with MadGraph at 



MCPM 



all. 
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Vertex 



Corresponding expression 




Table 3: Feynman rules for the most relevant vertices in the MCPM, part 1. 
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Vertex 



Corresponding expression 




Table 4: Feynman rules for the most relevant vertices in the MCPM, part 2. There are 
many more, including couplings of gauge bosons to Higgs bosons and Higgs-boson self- 
couplings. The arrows on the H ± lines denote the flow of negative charge. 
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Particles 



Names in MadGraph MCPM implementation 



Protons 






P 






Quarks 






u d c 


s t b 




Antiquarks 






u~ d~ 


c~ s~ t~ b~ 




Leptons 






e- m- 


tt- ve vm vt 




Antileptons 






e+ m+ 


tt+ ve~ vm~ 


vt~ 


Gauge bosons 






a z w+ w- g 




p', ti, h" 






rho hi h2 










h+ h- 






Table 5: The MCPM particles and their names in the MadGraph implementation. 


Section 


Parameter 


MadGraph 


Default 


Unit Strict MCPM 


Gauge 


l 

a 


aEWMl 




127.916 






G F 


Gf 




0.000011664 


GeV^ 2 




a s 


aS 




0.1184 






m z 


MZ 




91.1876 


GeV 


CKM matrix 


^WS 


lamWS 




0.2253 







Aws 


AWS 




0.808 






Pws 


rhoWS 




0.132 






Vws 


etaWS 




0.341 




Fermion masses m c 


MC 




1.42 


GeV 




m t 


MT 




172.9 


GeV 




m b 


MB 




4.67 


GeV 




m e 


Me 




0.00051100 


GeV 






MM 




0.10566 


GeV 




m T 


MTA 




1.7768 


GeV 


Higgs-boson masses m p i 


mrho 




125 


GeV 




m h > 


mhl 




200 


GeV 




mu" 


mh2 




150 


GeV 




m H ± 


mhc 




150 


GeV 


Table 6: The MCPM parameters 


and their 


names in the MadGraph implementation. Note 


that some parameters are zero in the MCPM (see 


last column), however they can be set to 


other values in 


this model. This allows the correct simulation of SM background processes 


(see section 3). 


The CKM matrix 


is given in the Wolfenstein parametrisation [31]. 
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